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In the Mott-insulator regime, two species of ultracold atoms in optical lattice can exhibit the 
low energy counterflow motion. We construct effective Hamiltonians for the three classes of the 
two-species insulators — fermion-fermion, boson-boson, and boson-fermion type — and reveal the con- 
ditions when the resulting groundstate supports super-counter-fluidity (SCF), with the alternative 
being phase-segregation. We emphasize a crucial role of breaking the isotopic symmetry between 
the species for realizing the SCF phase. 



Recent pioneering experiment by Greiner et al. |j| ( 
proposed in ref. on the superfluid-Mott-insulator 
transition in the optical lattice of ultracold atoms has 
opened up a new era of strongly correlated atomic sys- 
tems. Now it is possible to experimentally study a Mott 
phase of ultracold atoms. As shown in ref. ||] , the spinor 
(S=l) bosons in the lattice can exhibit very rich and ex- 
otic phase diagram. 

In the Mott phase of single specie of atoms, any low- 
energy transport is suppressed [Q. The situation may 
change dramatically, if the Mott phase includes al least 
two species. While the net number-of-atoms transport 
is still suppressed, the counterflow (when the currents of 
the two species are equal in absolute values and are in 
opposite directions), generally speaking, survives, and, 
at certain conditions, it can be non-dissipative (super- 
counterflow) . 

In this Letter, we discuss typical cases of counter- 
flow dynamics in two-species Mott phase. We consider 
three classes of the commensurate two-species systems: 
fermion-fermion (FF), boson-boson (BB), and boson- 
fermion (BF) mixtures (assuming that each of the species 
does not possess internal degrees of freedom). We con- 
fine our analysis to the strong coupling limit, when the 
effective counterflow Hamiltonian is readily obtained in 
the second order perturbation theory in the parameter 
t/U <C 1, where t is the hopping amplitude and U is the 
on-site repulsion. For the same-statistics (BB or FF) 
mixtures, the effective Hamiltonian can be written in 
terms of s = v jl iso-spin operators (the integer v is the 
filling factor). In the FF case, we have a s = 1/2 anti- 
ferromagnetic asymmetric effective Hamiltonian. In the 
BB case, the counterflow term is represented by asym- 
metric ferromagnetic interaction of iso-spins. 

Depending on particular means of breaking the iso- 
topic SU (2) symmetry, either the effective Hamiltonians 
are of the easy-axis type, which means phase segregation, 
or they are of the easy-plane type resulting in the SCF 
groundstate. In the BF case, it is possible to eliminate 



bosonic variables and obtain a spinless fermionic model 
with nearest-neighbor interaction. The sign of the inter- 
action depends on the parameters of the original Hamil- 
tonian, and how they deviate from the exact fermion- 
boson symmetry. In particular, the p-wave Cooper pair- 
ing is possible which results in the SCF. 

Our lowest-order effective Hamiltonians will involve 
only on-site and nearest-neighbor interactions. Hence, 
to derive them it is sufficient to consider a system con- 
sisting of two strong traps, 1 and 2, between which weak 
tunneling is allowed. Only the lowest one-particle state 
(p a for the species a =T>J. (we use pseudo-spin nota- 
tion) in each trap is taken into account. We introduce 
the creation and annihilation fcrmionic/bosonic opera- 
tors at , aba, respectively, where b = 1,2 labels the 
trap (site). Our original two-site Hamiltonian then reads 

#12 - H$ + V 13 , (1) 

Hlf = ^ ^2 Uaa ' :n b<T n b<7 r - , (2) 
b ctct' 

Vxi = - ML a 2 CT + H.c.) , (3) 

a 

where n blJ = al a a ba , U aa > = g aa > J (2x|¥v| 2 |</v| 2 , 9aa' is 
the interaction constant between species a and a', and 
:(. . .): denotes normal form of the product of the creation- 
annihilation operators. We consider the tunneling matrix 
element t„ as real and positive and depending on the 
component index. 

At commensurate total filling n b = n ba — v on 
each site, and in the limit \t a \ <C U aa >, single particle 
jumps change the total on-site populations, and, there- 
fore, require high energy. In contrast, exchanging two 
different particles does not require such an energy. In or- 
der to describe these processes the single particle jumps 
should be eliminated in the second order with respect 
to the ratio t/U. Thus, we choose in eq. ^ as 

zeroth order part, while Vn is the perturbation. The 
basis in our effective Hilbert space is given by the states 
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I off) = 1 7i if, nii; n2f) ^21), where ni = 712 = f ; the 
virtual excited states |ex) (where ni = n 2 ± 1) are 
excluded in the second order. The effective Hamiltonian 
V{ 2 can be represented in terms of the matrix elements 
of the original Hamiltonian as 

^ = -E^r|^p. ( 4 ) 

where a, (3 label |eff), and 7 denotes |ex); E's are 
the eigenenergies of H 1 - * 1 in eq. (||). 

FF case. This case is the easiest for calculation. The 
only possible filling factor is v — 1 and the only possible 
intermediate virtual state 7 is the state with two differ- 
ent fermions on one site. The only relevant interaction 
vertex is U^i = U. The Hamiltonian of the system is the 
standard Hubbard model j^] widely used for electrons in 
metal, the only specific (and significant !) feature being 
the dependence of the hopping amplitude on the "spin" 
index. The effective Hamiltonian is readily obtained from 
eq. (||), in terms of iso-spin operators 

S 6 = (1/2) £ olW^flto , (5) 

<7<j' 

where (cr) CT(T ' are the Pauli matrices. The operators 
Sh = (Sbxi Sby, Sbz), where the second index refers to 
the iso-spin components, obey standard spin commuta- 
tion relation at same site and commute with each other 
on different sites. The resulting effective two-site Hamil- 
tonian H' 12 — 2JS 1 S 2 + 2J' Si z S 2z extended to a lattice 
is 

Hff = [JSjSj + J'SizSjz] — BSi Z , (6) 

J = 2t 1 t L /U, J' = m~t L ) 2 /U . (7) 

where (. . .) stands for the pairs of nearest neighbors; 
B = jU| — fj,i, and [i a is the chemical potential of the o~- 
th component. If the two species are, e.g., just different 
isotopes of one and the same element, then the condi- 
tion ^ naturally occurs from the difference of the 
masses. 

The groundstate of this system (see, e.q., B) at 
B = is the spin-1/2 easy-axis anti-ferromagnet, given 
J > 0, J' > in eq. (0). This corresponds to the checker- 
board insulator (Neel state). If J' = and for infinites- 
imally small B 7 this state is unstable with respect to 
forming the canted phase, that is, the easy-plane state in 
which the staggered magnetization resides in the (x, y)- 
plane, with the ferromagnetic component ~ B oriented 
along the z-axis. This state maps into the superfluid 
bosonic system (see in, e.g., H and also [pi), which is 
equivalent to the SCF. For J' ^ 0, the canted state 
occurs by the first order transition at \B\ > B c , where 
B c m <J J'(2J + J'). In the canonical ensemble, when the 
numbers of atoms of both components, iVj and N±, are 



fixed, the SCF phase corresponds to a general case, while 
the checkerboard insulator occurs only at iVf = iVj_. As 
shown in ref. J8|, in canonical ensemble, the checkerboard 
insulator state can coexist with the superfluid, with the 
amount of the superfluid determined by iVj — iVj_ ^ 0. 

It is important to stress, that in our situation — in con- 
trast to the case of real easy-plane anti-ferromagnet — the 
staggered magnetization vector fundamentally cannot be 
attached to any direction in the easy plane due to the ex- 
act conservation of the atoms of each component. This 
guarantees the exact U(l) symmetry. As discussed later, 
the easy-plane U(l) symmetry can be broken only spon- 
taneously by forming the superfluid counter-flow vacuum. 

BB case. In this case we also introduce the site iso- 
spin operators according to eq. (|5|) with bosonic creation- 
annihilation operators. The expression for the effective 
Hamiltonian in a general case of non-equal interactions 
is rather cumbersome. It becomes much simpler in two 
important particular cases: (i) when U aa i = U + 5U aa i , 
with |5f/(ro-/| <C U, and (ii) when v = 1. We start with 
the former case in the limit \5U arT i\/U — ► 0. Then, the 
denominators in eq. (^J) are the same and equal to U . 
Hence, the intermediate virtual states 7 do not involve 
any projection operators, and the answer for the two- 
sites can be written as a bilinear form of s = v/2 iso-spin 
operators as H' 12 — — JSiS 2 — J 1 S\ Z S2 Z — Bvz{S\ z + 
S 2 z) ,B 12 = (u+ - tfj/U, where the expressions 
for J and J' are given by eqs. (|). When \8U aa .\ < U, 
the parameters J, J', and Bi 2 remain, to the first or- 
der in ISUcrer'l/U, the same. The main correction to the 
Hamiltonian is associated with the term H }° (ITJ) . Then, 
extending the two-sites Hamiltonian H' l2 + H^ 2 to the 
lattice and omitting trivial constant terms, we find 

Hbb = — J]] [JSjSj + J'Si Z Sj Z ] 

+ Y J [(D/2)(S l zf -BS lz ] , (8) 

i 

B = ^- H -2(u-l)(U n -U u )+pB 12 , (9) 
D = U u + U u - 2U n , (10) 

with p being the number of close neighbors on the lat- 
tice. Note that in contrast to the FF case, the minus sign 
now stands in front of J and J'. That is, if the D-term 
in eq.(^) is ignored, we arrive at the easy-axis ferromag- 
netic model, that implies a phase segregation. However, 
for large enough D in eqs.(P,^0|), the easy-plane ground 
state can be realized. The mean field condition for this 
(in the case B = ) is D > pj' 0. 

The effective Hamiltonian for the BB situation at v = 1 
can readily be found for arbitrary U aa i . The summation 
over 7 in eq.(||) in this case is very simple, because for any 
different states a and (3 in eq. (^), there is no more than 
one state 7 for which matrix elements differ from zero. In 
the diagonal terms a = (3, if two sites are both occupied 
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by the same iso-spin a bosons, the energy change due 
to jumping of either boson is U aa ; when the bosons have 
opposite spins, the energy change becomes U^. The fi- 
nal result acquires the form eq. (^), the parameters </, 
J', and B being 



J = 2t ] t i /U u , 

J 1 = -(t T + hf/Un + 2*?M T + 2tyu u , 

b 12 = 2(tyu n - tyu u ) . 



(11) 

(12) 
(13) 



At certain conditions, one obtains the easy-plane situa- 
tion J' < 0, that is the SCF. A solution for this can, in 
principle, be found exactly. We will, however, analyze it 
in a simplified situation when the system is almost SU (2) 
symmetric (like 87 Rb jn|), which means « U~n ~ 
fTj-j,, ti ~ t|, and will choose chemical potentials to have 
B = 0. We represent U n = U + U' + D/2, U u = 



U -U' + D/2, U- 



Tl 



u, t 



t + t' t 



t - t' 



with \U'\/U < 1, \D\/U < 1, \t'/t\ < 1. Then, 
J' < yields 



|t7t| < y/D/2U, D>0. 



(14) 



The easy-plane state is equivalent to the SCF. 
Now let us focus on why the easy-plane ground state of 
the models ([| ||) can support the superfluid counter- 
flow of the components. First, we note that, according 
to Holstein and Primakoff p[ , spin lattice is equivalent 
to the lattice bosons, with the spin commutation relation 
being essentially equivalent to the Bose commutation re- 
lation. Then, the site operators S„ = Snx + iS n y are 
proportional to the effective boson annihilation opera- 
tors b n pH , and the operator S nz — (v/2) — fyfi n . 
Accordingly, the formation of the easy-plane order pa- 
rameter 5+ = (S nx + iS ny ) = | *S" | exp(zc^) — (b n ) , 
where \S'\ ~ const ^ 0, is equivalent to the forma- 
tion of the Bose-Einstein condensate (BEC) (b n ) ^ of 
the Holstein-Primakoff boson b n . As long as the global 
invariance with respect to the phase ip of the bosons 
holds, which is insured by the conservation of the origi- 
nal species, and the bosons are interacting, the system is 
(counter-) superfluid. 

Let us obtain expression for the superfluid counter- 
current density I, which is the difference of the currents 
of the each specie. To be specific, we will discuss the 
BB case and keep in mind that the result for the FF 
case is, qualitatively, the same. In the long wave limit, 
the conservation law of the components takes the form 
S z + VI = 0, where ^(x) w S nz /£l is the difference of 
the densities of the species, with 51 standing for the unit 
cell volume and x being close to the site n. This con- 
tinuity equation follows from the Heiscnbcrg equation of 
motion iS nz = [S nz , Hbb] after performing the commu- 
tations and taking the long wave limit. Finally, replacing 
the operators by c-numbers, we find 



where d denotes the nearest neighbor distance; n sc f has 
a meaning of the effective superfluid density, with 
being the corresponding velocity. Thus, the easy- 
plane ground state of the models (f| ||) supports a super- 
counterflow of the two components. The general condi- 
tion for this is that the intrinsic symmetry (i.e., SU(2)) 
between the components is broken down to the U(l) 
group. 

BF case. Now we turn to the Bose-Fermi mixture on 
the lattice. The Hamiltonian (|l|) of two sites becomes 

H[°2 ] = J2 [(^o/ 2 ) Ub( - nb - !) + n b™b] , (16) 

b=l,2 

V\2 = — (ts a\a\ + t F c\c\ + H.c.) , (17) 

where at, ai, and el, q, stand for the site creation and 
annihilation operators of the bosons and the fermions, 
respectively; n& = a\ab, nib = c\cb] Uq > 0, U\ > 0. 
Below we will see that the SCF is possible, if (formal) 
symmetry between bosons and fermions (when Uq — 
U\, ts — tp) is broken. 

The effective two-site Hamiltonian can be written in 
terms of the fermionic operators only. The elimination 
of the bosonic operators is possible because the trun- 
cated Hilbert space is exhausted by different fermionic 
occupations, the bosonic occupations being unambigu- 
ously defined by the constraint + m& = v. Formally, 
one introduces new operators Cj = a\ci/y/v and ob- 
serves that, given our constraint, these satisfy standard 
fermionic algebra. [Note, that the vacuum for the new 
fermions corresponds to the Bose Mott insulator state 
with the filling factor v] The effective Hilbert space and 
the Hamiltonian are then automatically defined: 



Hbf — 



[-Kcjcj + (A/2) rniirij] - fi F V] , (18) 



(we have omitted the tilde in the fermionic operators). 
Here /ip = a denotes the resulting fermionic 

chemical potential, with fi being the difference of the 
bare fermion and boson chemical potentials (adding one 
fermion to the lattice implies simultaneous removal of one 
boson); the other constants are: 



k = 2 vtstF fU\ 



A = 2ti — 



1 



Ui 2U - Ui 



(19) 



M12 



t F + t 2 B v 2 t 2 B (v 2 -l) 2t%v(v + l) 



2U - U t 



U 



(21) 



I = UsdWip, n sc { 



Jpd 2 \S'\ 2 /(12fl), 



(15) 



and it is assumed that 2Uq — U\ > p2| . 

Note that in the case of the exact boson-fermion sym- 
metry (in the meaning defined above), A = and 
SCF is not possible, because the problem maps on the 
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ideal one-component fermion gas. Thus, the symmetry 
must be broken. From eqs. (|20|), it is seen that there 
is a region where A < 0. At, e.g., U\ — Uq, we find 
A = 2(t 2 F — t%)/U , and A < when t B > t F . Thus, we 
obtain one-component lattice fermions with the nearest 
neighbor attraction. This leads to the Cooper pairing 
(see, e.g., |13 ), which in our case implies the SCF. As 



our fermions are effectively spinless (spin-polarized), the 
pairing takes place in the p-channel, that is the super- 
conducting order parameter is characterized by broken 
inversion symmetry, and can exhibit broken time-reversal 
symmetry (see in Q ) . The particular structure depends 
on the lattice. 

Experimental realization and detection of the 
SCF. The essential condition for the SCF is that the sys- 
tem is in the MI regime with respect to the net-atomic 
transport. It is exactly the same condition formulated for 
the one component case . In the FF case at T = and 
\B\ > B c , the system will be in the SCF state. In the BB 
case, the three interaction constants U aa must be tuned 
in order to avoid the phase separation. This, however, 
can be achieved by creating slightly displaced lattice po- 
tentials |^] for the species in order to reduce the overlap 
of the states f, J, on each site, so that the value J7f| is 
reduced and the parameter D ( pLf ) becomes positive and 
large enough to insure the condition ( |l4"| ) for, e.g., 87 Rb. 

It is important to note that in the cases FF and BB, 
once T ^ the phase transition from the ordered (SCF) 
phase into the "paramagnetic" (i.e., normal with respect 
to the counter-flow) state may occur. A typical critical 
(Neel or Curie) temperature T c is given by the effective 
exchange constant as T c pj. Thus, the SCF requires 
low enough temperatures T < pt 2 /U. For typical pa- 
rameters employed in ref. j^], this reads as T < 10 _9 K. 
In the BF case, the SCF-normal counter-fluid transition 
temperature is given by the same energy scale. 

The SCF could be detected by observing non- 
dissipative exchange of the components through the opti- 
cal lattice. One of the possibilities consists of separating 
the components initially. Then, the species will be ex- 
changing their position in the oscillating manner. Rais- 
ing T above T c will result in the abrupt increase of 
the damping of these oscillations. 

The existence of the phase ip ([l5]) implies that 
the counter-vortex can be supported in the SCF phase. 
Then, the winding of the phase in the SCF, which has 
a boundary with the regular BEC phase of the BB com- 
ponents, will imprint this winding on the BEC phases, 
resulting in regular BEC vortices. Detecting these vor- 
tices will, then, signal a presence of the SCF-vortex. 

The order associated with the SCF is revealed in the 
fwo-particle density matrix. Its detection can be done 
by atomic scattering of fast atoms off the lattice |L5j . 
Specifically, the cross-section of the process, when in- 
coming fast particle of, e.g., sort f strikes the atoms 
in the lattice and "transforms" into the outcoming par- 



ticle |, is exactly given by the (iso) spin-spin correlators 
of the atoms in the lattice. This issue will be considered 
in greater detail elsewhere. 

Finally, let us discuss a formation of the SCF vortex 
in a rotating lattice. To be specific, we will discuss the 
BB case. The lattice rotation can induce such vortices, 
if the components have different masses M a . Indeed, 
the total mass current can be expressed as JS M ) = 
(Aff — M|)n sc fVy, where n sc f, ip define the counter- 
current (|l5|). External rotation of the lattice at some 
angular velocity uj changes the energy of the counter- 
vortex ~ n sc {L ln(L/d), where L stands for a typical 
system size, by ~ ±(Mf — Mi)n sc fL 3 uj. The formation 
of the vortex becomes energetically favorable at 



\oj\ >lo c k L- 2 hx{L/d)/\M^-Mi\. 



(22) 



Note that differences in masses can result in the phase 
separation induced by gravity/rotation. However, if the 
atomic magnetons of different species are different, the 
corresponding linear and oscillator fields can be compen- 
sated by adding non-uniform magnetic field. 

Summarizing, we have demonstrated that there is a 
large variety of strongly correlated groundstates in the 
two-component system of ultracold atoms in an opti- 
cal lattice, even if the system is in the Mott-insulator 
regime with respect to the net number-of-atoms trans- 
port. The most dramatic effect that may occur is the 
so-called super-counter-fluidity, when the system sup- 
ports a non-dissipative counterflow of the two compo- 
nents. In the strong-coupling limit, the effective Hamilto- 
nian for the super-counterfluids corresponds to: (i) easy- 
plane [canted) anti-ferromagnet (in the case of fermion- 
fermion mixture) , (ii) easy-plane ferromagnet (in the case 
of boson-boson mixture) , and (iii) p-wave superconductor 
(in the case of boson- fermion mixture). The experimen- 
tal conditions and means of detecting the super-counter- 
fluidity are outlined. 
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